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We study the coupled dynamics of the displacement fields in a one dimensional coupled-field 
model for drifting crystals, first proposed by R.Lahiri and S.Ramaswamy [Phys. Rev. Lett. 79, 1150 
(1997)]. We present some exact results for the steady state and the current in the lattice version of 
the model, for a special subspace in the parameter space, within the region where the model displays 
kinematic waves. We use these to construct the effective continuum equations corresponding to the 
lattice model. These equations decouple at the linear level in terms of the eigenmodes. We examine 
the long-time, large-distance properties of the correlation functions of the eigenmodes by using 
symmetry arguments, Monte Carlo simulations and self-consistent mode coupling methods. For 
most parameter values, the scaling exponents of the Kardar-Parisi-Zhang equation are obtained. 
However, for certain symmetry-determined values of the coupling constants the two eigenmodes, 
although nonlinearly coupled, are characterized by two distinct dynamic exponents. We discuss the 
possible application of the dynamic renormalization group in this context. 



I. INTRODUCTION 
A. Background 

Spatial and temporal correlations in spatially extended systems with a conservation law or a continuous invariance 
are widely observed to satisfy a scaling or homogeneity property. For instance, if the system is described by a single 
scalar field 4>(x,t), the correlation function C(x,t) = (<j)(0,0)(j)(x,t)} satisfies 

C(x,t) ^b 2x C{bx,b z t). (1) 

under rescaling of space by a factor b. Here z is the dynamic exponent while \ describes the spatial scaling of the 
field. Equation (1) holds in the rest frame of the ^-fluctuations, so that if the system has wavelike excitations, e.g., 
kinematic waves [Q of moving </>- fluctuations, it is necessary to perform a Galilean shift to co-move with the wave. 
The exponent z then describes the dissipation of the fluctuation wave, with a fluctuation of spatial extent Aa; having 
a lifetime proportional to (Ax) z . 

Now consider the scaling properties of systems with several coupled fields, say {<j) a , a — 1 to A}, whose dynamical 
evolution involves inter-field couplings both at the linear and nonlinear levels which result in propagating kinematic 
waves. At the linear level, the problem requires diagonalizing the matrix of couplings. The eigenvectors tp^ involve 
linear combinations of the <p a and represent modes that propagate as independent kinematic waves. The real and 
imaginary parts of the eigenvalues c M encode respectively the speeds and dampings of the corresponding kinematic 
waves, and in general, differ from one wave to another. By performing a Galilean shift with speed c M , one may move 
to the rest frame of mode ^; kinematic waves corresponding to other modes are not stationary in this frame, however, 
and these moving modes also contribute to the dissipation of mode /i as they are coupled nonlinearly to it. The 
correlation function C ti (x,t) = (^(0, 0)^^, t)} is expected to satisfy 

C^(x,t) « b 2x »C{bx,b z n). (2) 

where \n characterizes the spatial scaling of mode fj, and z^ is the corresponding dynamic exponent. 

The question arises: Is there a single common value z which characterizes the decay of all the modes fil When the 
answer is yes, as in fact it generally is, we say that the system obeys strong dynamic scaling. Considerable interest 
therefore attaches to exceptions to this general rule. Accordingly, one would like to characterize the conditions for 
the occurrence of weak dynamic scaling, when at least one z^ is different from the rest. A priori, there are two sets 
of circumstances when weak dynamic scaling may be expected: 

(i) If the transformation from <p a to ip^, which is designed to decouple modes at the linear level, actually succeeeds 
in decoupling them for the full nonlinear problem, then evidently each mode evolves autonomously and independent 
z^'s may arise. In fact, a complete decoupling at the nonlinear level does occur in the context of a reduced model 
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of magnetohydrodynamic turbulence |2|, and may well arise in other problems as well. In the MHD case, however, 
it turns out that both modes obey evolution equations with similar (autonomous) nonlinearites, so that a common 
value of z ensues. But this need not be the case for other problems. 

(ii) Consider a situation in which the evolution of a subset of the fields, say {u/3}, does not involve the others {u a } 
, while the evolution of the set {u a } does involve {up}. In this case, zp and z a may take on distinct values. Indeed, 
this is borne out by numerical studies of two-field dynamics || which show that weak scaling occurs if the evolution 
is autonomuous in one of the two fields, or very nearly so, in which case strong crossover effects may be expected. 

One of the principal results of this paper, which we present below, is a third possible scenario for weak dynamical 
scaling, where neither field is autonomous. 



B. Results 



In this paper we examine weak and strong dynamical scaling in a system with two coupled fields, which result in 
two coupled kinematic waves characterized by mode coordinates ipi and ip2 respectively. We work both with a lattice 
model involving two sets of spins, and with the corresponding continuum equations involving two coupled scalar fields. 
The analysis of the lattice model is facilitated by showing that along certain representative loci in parameter space, 
the steady state has product measure form. This allows the current to be found, and a continuum expansion to 
be made, with coefficients that explicitly involve the parameters and mean occupations of the lattice model. This 
enables us to make direct comparisons between the results of numerical simulations of the lattice model and analytical 
self-consistent calculations for the continuum equations. 

Our most interesting result is the identification of a third set of circumstances, beyond (i) and (ii) mentioned in 
Section I A, in which weak dynamic scaling results, despite each mode being nonlinearly coupled to the other. This 
involves symmetry properties of the kinematic waves: we find weak dynamic scaling if we choose model parameters 
so that the evolution equations are invariant under inversion of the second mode coordinate (ip2 - > — 1^2) but not 
under tpi — > —ip±. In that case our numerical simulations show that z\ = 3/2, while z 2 — 2 with multiplicative 
logarithmic corrections. We also study the problem within a self-consistent mode coupling calculation which shows 
that the different dynamic exponents arise in a rather interesting way: the linearized version of the problem has z = 2 
for both fields. The scattering of tjji by fluctuations in ip2 and vice versa gives rise to singular corrections to the 
diffusivity for ipi, leading to z\ — 3/2. The fluctuations, however, cause no singular correction to the relaxation of 
■02, leaving z 2 = 2. For most other parameter values, however, we find the more common strong dynamic scaling with 
z\ = z 2 — 3/2. 

The remainder of this paper is organized as follows. In Section |l|, we review briefly the continuum stochastic 
dynamical equations of Ref. Q. In section [ii] we present the lattice model, and show how the condition of pairwise 
balance can be used to find the exact steady state if the transition rates satisfy a certain relation. We characterize 
changes of the symmetry of the evolution equations as overall densities are varied, and report the results of extensive 
numerical simulations which show that the values of the dynamic exponents depend strongly on these symmetries. In 
particular, we present evidence for weak dynamic scaling when the two kinematic waves have different symmetries. In 
section ^ we describe analytical methods, primarily a one-loop self-consistent treatment of the continuum stochastic 
PDEs, for calculating the exponents in the weak dynamic scaling regime. We also outline a dynamic renormalization 
group procedure for this regime, discuss the difficulties that arise therein, and remind the reader how strong dynamic 



scaling is restored for generic values of the parameters in the model. We close in Section VI with a summary 



II. CONTINUUM STOCHASTIC PDES FOR DRIFTING CRYSTALS 

We review very briefly here the construction of our model equations of motion; details may be found in |^,^] . The 
physical system which inspired the initial work on the model was a lattice drifting through a dissipative medium. 
There are at least two examples of this: (i) steadily sedimenting colloidal crystals and (ii) a flux lattice driven, by the 
action of the Lorentz force of an imposed supercurrent, through a type II superconductor. If inertia is ignored, the 
equation of motion of the displacement field u(r,i) is of the form velocity — mobility x force, i.e., 

u= M(Vu)(DVVu + F) + C, (3) 

where the mobility tensor M is allowed to depend on the lattice distortion Vu, the tensor D represents elastic forces, F 
is the driving force, and ( a suitable noise source. Our results are for a highly simplified model with the same physics 
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as in (H). This model, constructed and studied in describes the coupled dynamics of two fields u x and u z (the 
displacements transverse to and along the drift direction respectively), as a function of one coordinate x transverse 
to the drift direction z. The equations of motion are 

u x = \izd x u z + -iid x u x d x ii z + Did x 2 u x + f x , (4) 

u z = M\d x u x + -f 2 {d x u x ) 2 + -iz{d x u z ) 2 + D 2 d x 2 u z + f z , (5) 

where f x and f z are zero-mean, Gaussian, spatiotemporally white noise sources. The equations arc invariant under 
the joint operations x — > —x, u x — > — u x . For A12 = A21 = 0, (0) and (Q) reduce to the Erta§-Kardar equations 
for drifting polymers with the larger symmetry x — > — a; (with or without iij — > The system can distinguish 

between up and down: there is no invariance under inversion of u z . The terms in (^) and (^) involving first spatial 
derivatives have the following interpretation: a tilt (d x u z ) produces a lateral drift (at a rate which depends on the 
density perturbations d x u x ), while the vertical speed depends both on compressions or dilations {d x u x ), as well as 
tilts (d x u z ). In this paper we shall consider only the case Ai 2 A 2 i > in which case the dispersion relation 

to = ±\Z\i 2 X 2 iq - iDq 2 . (6) 

holds for the linearized version of (Q) and (^) and predicts travelling waves at small wavenumber q. The linearly 
unstable case A12A21 < has been discussed extensively elsewhere ||. The mode coordinates corresponding to (||) 
are given by 

V>1,2 = VcUx ± \/fr"z (7) 

As discussed in ||, the long-time, large-lengthscale behaviour of the PDEs (f|) and <Q) are expected to be the same 
as those of a particular two-species Ising-Kawasaki model in which the jump rate of each species depends on the local 
density of the other. We turn next to this discrete model and its dynamics. 

III. THE LATTICE MODEL 

The lattice model is defined in terms of two sets of variables {a{\ and {r i _i} which reside on two interpenetrating 
sublattices with periodic boundary conditions; the former set occupies the integer sites and the latter the half- integer 
mid-bond locations of a one-dimensional lattice with L sites. Each cr, and r i— i is an Ising variable taking on values 
±1. They represent discrete versions of the density and tilt fields in the sedimentation problem: If ui is 1, there is a 
particle (+) at site i, and if o~\ = — 1, there is no particle (— ). The variable T i _i — 1 and — 1, implies two values / 
and \ of the local tilt respectively. A typical configuration of the full system is thus +\ — / — / + \ — / + / + / + \— . 

Both sets of variables are conserved, i.e. J2 a i an d J2 T i-± an d the associated densities p° = XXl + cr i)/2£ and 
p 2 = 2(1 + Tj_i)/2i are constant. We consider a r-dependent local field which guides the cr-current and vice versa. 
Thus, for instance, the Kawasaki exchange dynamics of the adjacent spins a% and ct^+i occurs at a rate which depends 
on r i+ i. The moves and the corresponding rates are depicted below: 
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The macroscopic behaviour of the model is determined by the relative values of the rates; a brief review of the 
phases and their characteristics appears in B. There are two distinct regimes separated by a nonequilibrium phase 
boundary as depicted in Fig. 1. The regime pi > p 2 , marked SPS in Fig. 1 was explored in detail in ||. In this 
phase, the system undergoes spontaneous phase separation of a particularly strong sort. Along the boundary pi = p 2 , 
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marked FDPO in Fig. 1, the system undergoes fluctuation-dominated phase ordering of a delicate sort, as discussed 
in 0. Finally, in the phase of the model with p\ < P2, marked KW in Fig. 1, there is no phase separation, and 
fluctuations are transported by kinematic waves. This is the regime of interest in this paper. 




FIG. 1. The phase diagram of the Lahiri-Ramaswamy model. The strongly phase separated phase (SPS) is separated from 
the kinematic wave phase (KW) of interest in this paper by the threshold line p\ = pi along which there is fluctuation-dominated 
phase ordering. The steady state can be found exactly along the line pijpx = rx/rz, by using the condition of pairwise balance. 

A. Exact results for the steady state 

The steady state can be found exactly provided that r\ = P2 and T2 — p\ (along the line marked Pairwise Balance 
in Fig. 1). To this end, let us choose new symbols to denote the values of the variables a and r: use 1 if a site or bond 
is occupied by a + or /, and use for — or \. Then the moves (l)-(8) listed above reduce to moves (a), (b), (a), (b) as 
follows: 

(1), (5) (a) : 100 — -> 001 p 2 (9) 
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The use of the new symbols 1 and explicitly brings out the fact that the dynamical moves on the two sublattices 
are alike for the choice of these special relations between the rates. This is crucial for further analysis. 

In this new representation, a configuration C is specified by the occupations of all sites. The time evolution of the 
probability P{C) of the occurrence of C is given by the master equation 

= W(C — ► C)P{C) - W(C — ► C')P{C) (13) 
c> 

where the transition rates W(C — > C) are given by p2 and p\ for changes of configurations involving moves (a) or 
(&), and (a) or (b) respectively. 

The dynamics preserves the sublattice densities p° and p^- Within the subspacc labelled by (pj,/^), one can see 
that the system is ergodic by noting that the dynamics induces a leftward drift of a pairs 00 and 11 with rate P2, and 
a rightward drift of such pairs with rate p\ . By successive application of moves (a), (b) in Eqs. @, © and their 
reverses, any configuration C in the subspace can be taken to a 'standard configuration' C Q in which all l's and all 
0's are clustered together. Since the lattice is periodic, the 11 and 00 pairs can be shifted to any other configuration 
C from C Q . Thus any configuration C can be taken to any configuration C via C a . 

In steady state the right hand side of ( |T3| ) must vanish. A sufficient condition for this is that fluxes balance in 
pairs, i.e. for every flux out of a configuration C to a configuration C , there should be an incoming flux from another 
uniquely determined configuration C into C. This is the condition of pairwise balance B: 
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W(C — ► C)P(C ) = W(C — ► C )P{C) (14) 

which is a generalization of the well known condition of detailed balance M . 

For our problem, C may be constructed as follows. Let us denote the configuration C symbolically as 
^m 1 Qm 2 im 3 Qm i ^m,5 , , . Qm t ^ w here there is a cluster of l's of size mi, followed by O's of length m,2, and so on, with a total 
of A; such clusters. Consider a transition to a configuration C = i m i() m2 ~ 1 l 2 01 m3 ~ 2 m4 l m5 • • • mk (a pair of 1 l's jump 
to the left, i.e. move (6)). One can always find a unique configuration e.g. C" = i™io m n m3 - 2 01 2 O m4 " 1 l m5 • • • mfc 
which gives rise to C (via move (6)). If the outgoing transition involves a rearrangement at the left edge of a cluster, 
the incoming transition involves a rearrangement at the right edge of the same cluster. Such an identification is 
possible also for transitions involving moves (a), (a) and (6), and ensures that W(C — ► C) = W(C — > C ). Thus 
Eq.(|l4"|) is satisfied provided the steady state probabilities obey 

P(C) = P(C") = constant. (15) 



This means that in steady state, every allowed configuration is equally likely. The constant appearing in Eq. (|1^) can 
be found on using the normalization condition P{C) = 1. If N-y = p\L and N% = p%L are the number of particles 
on the two sublattices, the total number of configurations in sector (p°, p%) is J\f = ( l Cn 1 ){ l Cn 2 )i where n Cm is the 
number of ways of choosing M out of a total of N objects, and hence P(C) = 1/W. 

In the thermodynamic limit L, N%, and N2 — > 00, with p°, p^ held constant, P(C) approaches the product measure 
form 

Y[p{o-Mr l+h ) = Rl Nl {l ~ pl) L - Nl fi N *(l ~ Pl) L - N * (16) 

i 

This form of the steady state holds also for a higher-dimensional generalization of the model involving rules (a) and 
(b) and their reverses along the sites and bonds in the d directions of a simple cubic lattice. 

The product measure weight in the steady state implies that correlation functions on different sites decouple. This 
then allows the current of a particles 

,,, . (K _ ,,)< <' +•»>(' -^fr-W , + ta -m^I^I^H) (,„ 

to be found explicitly, 

Ji = (p 2 -pi)pi(l-pi)(l-2p 2 ) (18) 

The first term in Eq. (0) comes from a particle hopping between sites i and i + 1 in the absence of a particle at site 
i + I while the second is for hopping in the presence of a particle on the site in between. A similar expression holds 
for the r-current J2, 

Ji = (p2-pi)pa(l-flj)(l-2pi). (19) 

The product measure form also allows us to find the roughness exponent of an associated height model, where 
the height fields associated with {at} and {rj} are respectively tin — J2k=i( a k ~~ ( a k)) an d = X)l-=i( T fe — ( T fe))- 
Fluctuations of the height field are characterized by the root- mean-square height difference G\ (r) = */ ((hu+ r — ha) 2 }, 
with G2(r) defined similarly in terms of {/i2i}- Using the fact that (eroCfc) = (&a){&k) for k ^ and 1 for k = 0, we 
find Gi(r) = G 2 (r) - r 1 ' 2 . Thus the roughness exponent \ defined by the growth of the root-mean-squared height 
fluctuations is 

X - 1/2 (20) 

for both height fields. Evidently, this value will also characterize fluctuations of linear combinations of the height 
fields h\ and ti2, which arise when we deal with mode coordinate fields in the next section. 
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FIG. 2. The decay of correlation functions away from the line of pairwise balance is shown for Cn (filled circle), C22 (filled 
diamond) and C'12 (filled triangle). When pairwise balance does hold, the correlation function vanishes (empty triangles). 

We close this section with numerical results for correlation functions away from the pairwise balance (PB) locus. 
We set r 2 = pi = 0, and investigate what happens if we move away from the PB locus r\ — P2- We studied the spatial 
correlation function C(r) = {/Ji<7i +r ) — (cr^) (<x; +r ) by Monte-Carlo simulation, and show our results for different values 
of n and P2 in Fig. |[ As expected, for the PB case n = P2 = 1 (empty triangles in Fig. |^), the correlation length 
is zero. Away from PB, with n = 1/3 and P2 = 1, the three curves for Cn, C22 and Ci 2 representing respectively 
the intra-sublattice correlation functions for sublattices 1 and 2 and the inter-sublattice correlation function all decay 
differently, but with a finite correlation length. (Fig. ^). This shows that although product measure does not hold 
away from the pairwise balance locus r\ — p 2 , correlations are short-ranged so that the behaviour on large length and 
time scales is expected to be similar to that in the pairwise balance case. 

B. Continuum equations and kinematic waves 



The expressions ( |18| ) and (|19|) for the current help us to construct approximate continuum equations for the evolution 
of the density fluctuation fields. The starting point is the pair of continuity equations 



dp m 
dt 



d_ 

dx 



Jm{pi,P2) 



1,2 



(21) 



where p\(x,t) and p 2 (x,i) are coarse-grained densities at a mesoscopic scale and Ji(x,t) and J 2 (x,i) are the corre- 
sponding currents. Each of these currents is written as the sum of three parts, 



J n 



rsystematic _r_ jdif fusive _i_ 
J m ' J m ' 'In 



(22) 



The systematic contributions j^v stematlc a t the mesoscopic scale are assumed to be given by the expressions ( |i~8| ) and 
( |l9| ) for the currents in an infinite system. The diffusive part J^ff usme arises from local density inhomogeneities and 
is taken to be —D m ^r. Finally r\ m is a noise term added to mimic stochasticity at the mesoscopic level; we consider 
uncorrelated white noise with (77m) — and (rj m (x, t)rj m (x' , t')) — T8{x — x')S(t — t 1 ). These continuum equations 
have the same symmetries as the lattice model and hence would be expected to exhibit the same behaviour on large 
length and time scales. 

Writing p\ — p\ + p\, and P2 = P2 + Pi (where p°'s are fixed average densities and p's are fluctuations) and using 
Eqs. (0), (|l9|), ^l| and [2^, one can write down coupled equations governing the evolution of p's. We write these 
in terms of the height functions hi = J p\{x ,t)dx and = / Pi{x ,i)dx , the continuum analogs of the discrete 
functions hu and /i 2 i defined in the previous section. We find 
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(23) 



and 



dh 2 , 



2p° x ){l-2p° 2 )^ 



dh 

2r ^(l-^-^+fl., 



dx 2 



(1-2P!)(^) 



,<%2 ^ 2 ,dh% 



(24) 



where r' = (fx — f2)> 

Let us define a = r {l-2p° 1 )(\-2p 2 \), b = r'p%(l-pf), c = r p°{l - p°), Kl = r'(l-2p°), and k 2 = r'(l-2p°). It is 
apparent that by taking linear combinations one can construct eigenmode fields h± = \fch\ ± \fbh 2 which decouple at 
the linear level. These fields describe wave-like modes JlJ travelling with speeds c± = —a±2y/bc. The time evolutions 
of these fields h± are governed by 



dh^ 
dt 



dh- 
~dt 



dh 



= c+ 
1 



+ 



9a; 

H 2 



dh. 



dx 



dx 2 

- D ff2h - 



3 , K 2 

1 



+ \2 



k 2 Ki . dh+ dh 



Vb 



dx 2 



1 , k 2 k% ..dh 
2 [ 7 = c~7b ){ ~dx 



( dh+ 
2^/bc dx ' 



,dh^ 
1 dx ' 



+ \2 



dx 



dx dx 
+ V+(x,t); 



2 1 ^ V6 ^ j + V 



+ \2 



1 ,dh 



2Vbc { dx ] 



(— ) 2 -( 

y dx' K 



<9.r 



<9x <9x 
+ r)-(x,t). 



(25) 



The new noise terms rj± — \fcr]i ± Vbrj 2 , are also delta-correlated. We have assumed D\ = D 2 = D, though this may 
not be preserved in the effective long wavelength equations. The fields h + and h- are coupled at the nonlinear level, 
so that each wave influences the dissipation of fluctuations of the other. We consider the dissipation properties of the 
waves in the next section for different sublattice filling fractions pi and p 2 of the sublattices. 

IV. DISSIPATION OF THE WAVES AND DYNAMICAL EXPONENTS 



If it happens that some of the coefficients of (g5|) vanish for certain choices of densities p° and p 2 , the evolution 
equations have special symmetries and this can have important implications for the long-time dynamics. As discussed 
below, there are three different symmetries which arise in the coupled field problem, each corresponding to a different 
set of dynamical exponents. The dynamical exponents associated with the wave modes, may differ. We have considered 
three special pairs of densities (p°,p 2 ) corresponding to three different symmetries. 

A. Symmetries of the equations: 



To facilitate subsequent discussions let us first consider the case of a single field h, and list four different symmetries 
]To| for its evolution. 

(a) RI symmetry: Invariance under up-down Reflection (i?) symmetry h — > —h and under Inversion (J) of space 
x —>■ —x. 

(b) RI symmetry: Invariance under h — > —h and not under x — > —x. 

(c) RI symmetry: Invariance under x — > — x and not under h — > —h. 

(d) RI symmetry: Invariance neither under x — > —x, nor under h — > —h. 

Since ^ is odd an equation of motion which contains only terms odd in h will be said to have R inversion symmetry, 
any term that is even in h will be said to break R symmtery. Accordingly: A term like ^7 obeys RI symmetry. 
Terms like and (§|) 3 obey RI symmetry. The RI symmetry is respected by the term (§|) 2 , while a term like 
(||) added to it breaks that and gives rise to RI symmetry. 

To illustrate the occurrence of different types of symmetries in our coupled-field problem, we consider three special 
pairs of densities (p°,p 2 )- 

(I) For pi — p 2 — h, (p5|) reduces to a pair of coupled equations, with linear and first and second derivative terms 
and cubic gradient nonlinearities: 
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(26) 



These equations describe two kinematic waves moving with speed c+ = r /2 and c_ = — r jl. The nonlinear couplings 
imply that each wave influences the evolution of the other. In order to study the dissipation of say the +mode, it is 
essential to move to the frame which co-moves with it. This is accomplished by a Galilean shift x — > x + c + t, and 
t — > t. In this frame, the evolution equations become: 
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(27) 



Evidently, in this frame, the —mode has a speed c_ — c+ 



x, because of the linear and cubic nonlinear terms. The RI symmetry 



The pair of Eqs. (|27|) are invariant under 



dx 



h+ — > — h+, h- — > — h- but not x - 
holds in the rest frame of h + mode. 

Similarly the dissipation of the —mode can be studied by going to a frame which co-moves with the —mode. It is 
easily seen that RI symmetry holds in this frame as well. 

Let us recall what happens when RI symmetry holds in the case of a single field h. In the Edwards- Wilkinson (EW) 
equation 11 1, the presence of an additional cubic term like (§§) 3 reduces the RI symmetry to RI symmetry. About 

the linear fixed point (with \ = 1/2), the cubic term has the same naive scaling dimension as 

( fpr) — both scale as 

6~ 3 / 2 when x — > bx and h — > b x h. Such a marginal cubic term is known |12Wlq] to introduce logarithmic factors in the 
behaviour of the height-height correlation functions. Using mode-coupli ng p2[ and dynamic al renormalization group 
p3[ calculations, it was found that the correlation function F(t) = yj ([h(x, t) — h(x, 0)] 2 ) grows as t l ^ A {log{t)) 1 ^ . 
Recalling that F(t) ~ t 13 with (3 — x/ z , we see that despite the lack of I symmetry, the critical exponents (3 and z do 
not change from their EW values 1/4 and 2, respectively. 

In our case with two coupled fields, the cubic gradient terms again have the same naive scaling dimension as the 
linear second order term, but they are more complicated than just (§§) 3 - We might guess nevertheless due to the 
symmetry that each of the + and —modes have z = 2. We will present numerical evidence in section 5.2, which 
confirms this and shows that there are similar multiplicative logarithmic factors. 

(II) For pi = 1/2 and p% ^ 1/2, the equations for the mode fields reduce to: 



dh+ 
~dT 



dh 



= 2Vb~c—~ 
dx 

1 ki ,dh 



■dh+ d 2 h 



2yfc y dx ' 
dh- 



<+ 

dx 2 
1 



3 n 2 ,dh + 2 
2 ^fc ( ~dx~ ) 



K2 dh + dh- 



c dx dx 

' ( dh + ]2 _ ( 9h-, 2 
[ dx ' [ dx ' 



>h- „ r— 



dx 



- v 



1 k 2 dh+ 2 i 

2^ [ dx> 2Vb? [ 



2Vbc dx 

d 2 h- 3 k 2 ,dh_. 9 K2 dh + dh 
1 -( ) H - — — — 

dx 2 2 y/c dx y/c dx da 

I ,dh- s \ ^dh + , 9 t dh- 



r)+(x,t); 



dx 



-) 



yfc dx dx 

ri-(x,t). 



dx 



(28) 



Going to either of the frames in which the -(-mode or the — mode are at rest, we see that with the cubic nonlinearities 
the Rl symmetry applies for each of the fields; R symmetry is broken by by quadratic nonlinear terms, and I is 
broken because of linear first order and cubic terms. The most relevant terms at the linear fixed point are the 
quadratic nonlinear terms. Thus we would expect that these terms would govern the dissipation and give rise to the 
Kardar-Parisi-Zhang (KPZ) value z — 3/2 jl6| for both the modes. 

(Ill) For p\ = 7^ 1/2, we have k\ = K2 and c = b, and the following pair of equations hold: 



dh+ 
dt 



dh+ 

dx 
K2 ,dh- 
~Tc^~dx~ 



D 



d 2 h^ 
dx 2 



3 K 2 ^dh+ 2 
Jc dx 



( dh+ 



2\/bc dx 



,dh 



dh- 



dx K *~ 1 



dx 



■r) + (x,t); 
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dh_ <9/i_ ^d 2 h_ Kodh+dh 
= c-—+D-— + 2 ' 



dt dx dx 2 \fc dx dx 

+ r)-{x,t). (29) 



i at± 



2Vbc dx 



( dh + ]2 _ ( dh --,2 
1 dx ' 1 dx ' 



Here, an interesting situation arises. In the co-moving frame of the —mode, the pair of equations are invariant under 
h_ — ► — h- and h + — > h + but not under x — * — x. Thus the /i_ field has i?/ symmetry, while the moving /i + field 
has RI symmetry. The equations are invariant under /i_ — ► — /i_ in any frame. The invariance under x — ► — x is 
broken (a) by the linear and (b) by the trilinear terms. The effect of (a) can be shifted away by comoving; the 
(b) terms are expected to provide logarithmic multiplicative corrections. The same symmetries hold in the rest frame 
of the + wave. Based on these observations, we expect z — 2 for the —mode (perhaps with multiplicative logarithmic 
corrections), and z = 3/2 for the +modc. This is precisely the weak dynamical scaling discussed in the Introduction. 

The quadratic nonlinear terms in Eq. ( p9[ ) are exactly like those obtained in PJT^|- Similarly the nonlinearities 
in ( p8| ) are of the form obtained in |Q. The crucial difference is that we have additional linear gradient couplings, 
which bring in the possibility of observing weak dynamical scaling in our coupled-field system. We now turn to the 
next step of checking numerically and analytically our symmetry-based expectations. 

B. Growth exponents from Monte-Carlo simulation: 

For numerical simulation, we used a definition of height hi(t) which differs slightly from hi(t) discussed so far. 
Instead of defining heights as density integrated over space with respect to a fixed site, we define hu and h 2 i as 
integrated densities but with respect to the first particle, which is itself moving. Such a definition was used earlier in 
p9| , and was found to markedly reduce the fluctuation in the height-height correlation function. If the particles on a 
particular sublattice are labelled 1, • • • , Np, then the height hi(t) = rii(t) — (ni(t)) where rii(t) is the tag number of the 
particle at site i and the subtracted part has a linear time dependence. If the site i is empty, hi is determined from the 
tags of the closest particles on either side by a lever rule. If i (t) (whose average value = vt with v — particle speed, 
is the location of the particle 1) one has hi(t) = Ylk=i (t) P k ~ P°^o(t), using the fact that rn{t) — J2k=i (t)(P k + P°)- 
In our problem, the height-height correlation (hiA(t) — ^^(0)) grows as —VAP°t, where va = (1 — — 2^>g) • Th° 
corresponding continuum equation for h\(x, t) has two more additional terms compared to (p3[), i.e. a constant — vaP° 
which can be removed by an appropriate shift, and a noise term £ = VAtpi(x, t). The scaling dimension of C is lower 
(£ — > b~ z () than the noise rji (rji — ► b^ z ^ 2 ^ 1 ^ 2 r)i), and hence is less relevant. 

We are interested in knowing the dissipation property of the ± wave modes, so we numerically compute the height- 
height correlation function of the variable h± in the co-moving frame of the relevant mode, i.e. by a Galilean shift 
with speed c±. The new fields h± are the counterparts of h± defined earlier. We monitor the correlation functions 

S 2 z {t) = { h±(x- c±t,t)-h±(x,Q) + J£t + p°±c±t 2 ) (30) 

where the average (• • •) is both over space and the ensemble of configurations in the steady state, while J± = 

^fcvAP\ ± VbvsP2 an d p± — \fcp\ ± Vbp2- We expect S±(t) to grow as t 13 with f3 = xl z - 

We used lattice sizes L = 54000, while the number of history averages for each case is mentioned in the figure 
captions. There are two curves in each of the figures ||, |^, and ^| below: for example in Fig. |^(b) the lower one 
represents S 1 !, and the upper one represents the height-height correlation function of the moving —mode in the rest 

frame of the +wave, namely St = ( h_ (x — c + t, t) — h_ (x, 0) + J°_ t + p _c + t ) . Similarly S\ in Fig. |^(a) is defined 

in the rest frame of the —mode. The correlation functions S 2, s increase linearly with time as they sense the effect 
of a moving kinematic wave. In a discrete lattice simulation, x in Eq. (|3^) gets replaced by discrete integers i. The 
kinematic wave speeds are nevertheless real fractions, so height-height correlation functions S 2 and S 2 have oscillations 
of the period l/c±; these oscillations are noticeable for low t, but their relative contribution dies down for larger t. 
(I) We first discuss the height-height correlations for pi = 0.5 and p^ = 0.5 (Fig. ||) corresponding to RI symmetry. 
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FIG. 3. Height correlation functions for p° — p% = 0.5 (RI symmetry) averaged over 6 histories, (a) In the rest frame of the 
— mode, S__ grows as ~ t 1 ' 8 implying /3 = 1/4, while ~ t, (b) In the rest frame of the + mode, grows as ~ t 1/2 , while 



The kinematic waves go around the whole system with a time period t a = \L/c±\. With L = 54000 and c± = 
±r'/2 = ±1/4 we have t a — 216000 as r' — 1/2 in our simulations. The curve for the mode which is moving in the 
rest frame of the other should dip to a minimum at time t a , and that is why it shows a flattening at times around 
t = t Q /2 = 108000 when the moving wave has travelled halfway around. 

Figure ^ gives strong evidence for (3 = 1/4 for both modes in their respective rest frames, when p° — p% — 1/2. 
Thus z = x/0 is equal to 2 for both modes. More careful examination of the data reveals that there are actually 
multiplicative logarithmic corrections to the leading power law behaviour. The data is consistent with a growth 
S(t) <~ t 1 ' 4 (/o5(t)) 1 ' 8 , as may arise from cubic nonlincarities as discussed above. Evidence of this is seen in 
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FIG. 4. Effect of scaling the rest-frame height correlation functions of Fig. | by a factor R(t) = t 1 ^ 2 (upper curves) and 
R(t) = t 1/2 {log{t)) 1/A (lower curves). The flattening of the latter confirms the presence of the multiplicative logarithms in both 
Si and S%. 

Figs. ^(a),(b) which indicate that the data of Figs. ||(a),(b) have factors (logit)) 1 ^ multiplying i 1 / 2 . As discussed 
in section 5.1, the latter are probably due to cubic nonlinearities. 

(II) For p° = 0.5 and p% = 1/3, the system has RI symmetry. Figure || shows that both and S 2 . grow as ~ i 2 / 3 
implying that (3 — 1/3 and z — 3/2. Thus the dissipations of both the waves are KPZ-like. 
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FIG. 5. Height correlation functions for pi = 1/2, P2 = 1/3 (RI symmetry) averaged over 4 histories, (a) S_ grows as ~ t 
implying (3 = 1/3, while ~ t, (b) too grows as ~ t 2 ^ 3 and S 2 . ~ t. 

(Ill) Finally Fig. ^ shows the height-height correlation functions for p° = 1/3 and p§ = 1/3 (RI symmetry for one 
mode, RI symmetry for the other). 
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FIG. 6. Height correlation function for p° = p% = 1/3 (RI symmetry for one mode, RI symmetry for the other) averaged 
over 12 histories, (a) St grows as ~ t 1/2 and S% ~ t. (b) S 2 . grows as ~ i 2/3 and 5*i ~ t. 

We see that grows as ~ t 1 !" 1 with an indication of multiplicative logarithmic corrections, indicating 0=1/4 
while 5^ grows as ~ t 2 / 3, implying f3 = 1/3. Recalling that \ = 1/2 for both modes, we have z = 2 for the — mode 
and z = 3/2 for the + mode. There is a logarithmic factor (log(t)) 1 ' 4 ', multiplying t x l 2 for the —wave as is apparent 
from the flattening of the curve on dividing S?. by t 1/ ' 2 (Zog(t)) 1/ ' 4 (Fig. [7]). 
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FIG. 7. Effect of rescaling Si in Fig. |by R(t) = t 1/2 (upper curve) and R(t) = t 1/2 (Zog(t)) 1/4 (lower curve). The flattening 
of the latter confirms the presence of the logarithmic factor in S 2 . 

It is remarkable that although the waves are coupled to each other, two different dynamical exponents arise in the 
same system, in conformity with our expectations on symmetry grounds. This is the first instance we know of in 
which such a property arises in a fully coupled system in which neither field evolves autonomously. 
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V. ANALYTICAL DEMONSTRATION OF WEAK DYNAMIC SCALING 



We argued in Sec. [0] in the context of the model equations 29, and showed by using Monte Carlo simulations, that 
there are two different dynamical exponents 3/2 and 2 respectively, for the two eigenmodes h+ and h- in the model. 
In this section we use self-consistent mode-coupling and renormalization group methods to study the large-distance, 
long-time properties of the correlation functions: We calculate the roughness exponents \i,2 and the dynamical 
exponents Zx$ of the fields Our analytical results agree with the previous numerical results. For appropriately 

chosen values for the densities of the particles in the lattice model [see above eqn. (p9|) 1, the continuum equations in 
the center-of-mass frame take the form 

h+ - ^-d x h+ + ^-{d x h + ) 2 + y(cU0 2 = v + d xx h + + /+, 

h- + ~Y 9 xh- + \ 3 (d x h + )(d x h_) = v_d xx h_ + /_, (31) 

where B a /2 = c + = — c_, Ai/2 = Zk^/^/c, A2/2 = —K2/2, A3 = 2nil^Jc- We denote the dissipation coefficients by V\ 
and vi (we use two separate symbols in anticipation of weak dynamic scaling) . We also ignore the cubic nonlinearities 
for simplicity. The waves can be removed from either Eqs.(|3l"|) separately by comoving with the left and right going 
waves through the Galilean shifts: x — > x + ^f-t, h + — » h + and x — » x — h- — > h- respectively. In the comoving 
frame of h + Eqs. ( |3~l| ) become 

K + ^j{d x h + ) 2 + y (d x h_) 2 = v+d xx h+ + /+, 
h- + B d x h- + \ 3 (d x h+)(d x h-) = v^d xx h^ + /_, (32) 
whereas in the frame of h- Eqs.(|3l|) reduce to 

h+ - B d x h h + ^-{d x h + ) 2 + y (d^-) 2 = v+d xx h + + /+, 

h- + \ 3 (d x h + )(d x h-) = v-d xx h- + /_, (33) 
The bare response functions of the two fields h + and /i_ are given by 

G+ (k, u) = . 1 ; G- (k, w) = \ (34) 

in the right moving frame and 

G$%u)= ■ , • n 1 1 , 2 ; G-(k,u)= . 1 (35) 
in the left moving frame. The noise correlations in both the frames are given by 

(/ +) -(0,0)/+,-(M)) = 2D+,-5(x)6(t). (36) 

Note that there is no frame in which the drift terms in both equations vanish. Noise correlations however do not 
change with change in reference frames as they are delta-correlated in time. 

In this section we analytically calculate the dynamic and roughness exponents \ii z ii i = +, — for the fields hi, i = 
+,— respectively, defined by ([hi(x,t) — /ii(0,0)] 2 ) = x Xi g(x Zi /t). There are two limits in which their behaviour 
is well-understood analytically. First, in the absence of the kinematic wave (B a ) terms, Eqs. ( |3l| ) reduce to the 
EK H equations. Although the complete phase diagram of the latter is not known, they do have a locally stable 
renormalization-group fixed point belonging to the universality class of the KPZ equation [Hf l, with x+ — X- = 1/2 
and z + = Z- = 3/2. Secondly, in the absence of the nonlinear terms coupling h + and the kinematic wave terms can 
be removed separately in each equation by opposite Galilean transformations, yielding scaling properties independent 
of the wavespeed. In the previous section our Monte Carlo results show that even when there is a nonlinear coupling, 
for particular densities, there is weak dynamic scaling: z + = 3/2, z_ = 2. Below, we offer an understanding of weak 
dynamic scaling in an analytical framework. 

The possible occurence of weak dynamical scaling in this problem, makes a renormalization-group treatment difficult: 
one cannot rescale time in two different ways for the two fields. We therefore adopt a self-consistent integral-equation 
approach, for which such weak scaling presents no difficulties. We then discuss how to circumvent the problems in 
applying an RG treatment. 
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A. Self-consistent mode coupling calculation 



1. KPZ equation in a moving frame 



Before embarking on a calculation for our model, it is instructive to look at the simpler case of a growing KPZ 
surface in a moving frame. The KPZ equation will be supplemented by a linear first order gradient term: 



dh 



dh 
dx 



iA7 2 /i + f], 



(37) 



with (r)(x, t)r)(0, 0)) = 2D5(x)5(t). We ask: Does the wave affect the exponents? We know \ = 1/2 for the KPZ 
equation in 1 dimension. Since x gives the static probability distribution, it is independent of reference frames. Due 
to the exponent identity \ + z — 2 JlQ] , it follows that z = 3/2. Thus in this particular case exponents are not affected 
by waves. The question is, can we see this in a self-consistent or renormalization-group calculation? Recall that the 
dynamic exponent is given by the width of the peak of the dynamic structure factor. The location of the peak, as a 
result of the wave induced by our transformation to a frame moving with speed c, is at lu = ck. The response and the 
correlation functions in that frame are thus given by G (k,iu) = iuJ _ i( ! k+uk i and C (k,to) — ^-cky^+^k 1 res P ec tively. 
The one-loop integral that produces a singular correction to v is 



dqdQ 



2Dq(k - q)kq 



p - ck) 2 + v 2 q 4 ][i(uj -Q)- ic(k -q) + v{k - q) 2 



dq 



2Dq(k-q)kq 
vq 2 \uq 2 + v(k — q) 2 



(38) 



which is same as that one obtains in a standard calculation for the KPZ equation in the rest frame. A similar 
expression holds for the correction to the correlation function. In both these integrals, the wave can be shifted away 
trivially. Thus we obtain z = 3/2 and x — 1/2, in agreement with our expectations. With this background we now 
present the calculations for the exponents for our model. Interestingly, we will find that the effects of the waves cannot 
be trivially shifted away as they could in the simple example discussed above. 



B. Self-consistent calculation for the coupled model 

The mode coupling approach to solving equations such as (|3l]) consists in obtaining diagrammatic perturbation 
expansions for the renormalized propagator G and correlation functions C and rcsumming these in such a way that 
all the internal lines are renormalized correlation functions or propagators. This provides an exact solution for 
asymptotically small wave number q when the vertex corrections vanish for q — > or due to some fundamental 
symmetry in the problem. In that case the problem is reduced to solving nonlinear integral equations for G and C 
whose order is the same as that of the nonlinearity in the modified Langevin equation. 

In a one-loop self-consistent modecoupling theory, one writes down one-loop intergral equations for the response 
and correlation functions. The basic assumption is that there are one-loop corrections which diverge in the infrared 
limit. In a theory where there is no vertex renormalisation due to some Ward-Takahashi identity (arising from some 
continuous symmetry of the equation of motion) these equations are exact, because any higher loop corrections can be 
incorporated within one-loop dressed response or correlation functions. However, when there is no such symmetry of 
the system to prevent vertex renormalisation the above assertion is not true. In such a situation, one also has to write 
down a one-loop self-consistent equation for the vertex, which now is to be solved simultaneously with the one-loop 
equations for response and correlation functions. An example of this has been worked out in |po| . Here we work with 
Eqs.©. There are no diagramatic corrections to the Ai vertex at the one loop level. However, there are no such 
conditions on A2 and A3. In our particular problem, however we show using bare response and correlation functions, 
that one loop vertex corrections for A2 and A3 are all infrared finite. Hence we ignore them in our calculations. In 
Figs. |^ and [)]we show the one loop diagramatic corrections to A2 and A3 respectively. We first work in the comoving 
frame of h + : In that frame C+(A;, lu) has a peak at lu = B Q k and C-(k, lo) is peaked at lu = 0. 
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FIG. 8. One- loop diagrammatic corrections to A2; A line indicates a response function and a line with a small filled circle 
refers to a correlation function. 
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FIG. 9. One- loop diagrammatic corrections to A3; A line indicates a response function and a line with a small filled circle 
refers to a correlation function. 

Let us look at one of the diagrams very carefully (first diagram in Fig.^) : 

J iM ^/o f dndqD + q 2 q2 2zfcA 2 A| , 1 1 / 

I = ik\ 2 \ 3 2 — 5T = — — + — ) dqD+. (39) 

7 [(U - B q) z + vtq^yW + v+q*] B A Q v+ V- J 

This has no infrared divergence. Similarly, all other diagrams in Figs.(|^) and (|9|) are also finite in the infrared limit. 
We ignore all these finite corrections for A2 and A3, i.e., we ignore vertex corrections in our self-consistent analysis. 
We again justify this a posteriori, by showing that the vertex corrections remain finite in the self-consistent theory. 
The time-displaced correlation functions of the fields h + and h- arc 

C+(r,t) = (h + (0,0)h + (r,t)), C_(r,t) = (fe_(0,0)/i_(r,t)). (40) 

The scaling forms for the correlation and response functions as a function of wavenumber k and frequency lj are 

C + (fc,^) = fc^- 2 *+^+/ + (fc z +M, C-(k,L,) = k- 1 - 2 *—*- f-(k*~/u), (41) 

and 

G+(k,u) = k- z +g + (k z +/uj), G-(k,w) = k~ z - g_{k z - j uj). (42) 

Here, z + are z_ are the dynamic exponents and \+ an d X- are the roughness exponents of the fields h+ and h—. 
Notice that we have allowed the existence of two different dynamic exponents. Since the fields are decoupled in the 
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linearized theory, there is no cross propagator. The following one-loop diagrams contribute to the respective self- 
energies £_|_(fc,a;) and E_(fc, oS) (defined by G" 1 = G^ 1 — Ej, where G^o is the bare propagator for hi, « = +,—) of 
the fields h + and h-. 





FIG. 10. One-loop diagrams contributing to E_(fc,u). A line refers to a response function and a line with a small filled circle 
refers to a correlation function. 



Notice that G^ik^uj = 0) = iB k + v + k 2 , GZq = v-k 2 . A self-consistent calculation is required if one encounters 
infrared divergences in the bare perturbation theory. The following one loop diagrams in Figs. 10 and |ll] contribute 
to E_(fc,w) and G_(fc, u>) respectively. It is easy to see that none of these diagrams diverge in a bare perturbation 
theory. The first diagram in Fig. 10 has the form 



E_(fc, ui) ~ A2A3A: / dq 



D_q 2 (k-q) 



v + q 2 [v + q 2 + v-{k - q) 2 - iB (k - q)] 



A2A3A: 



A2A3A: 



D^q 2 (k - q)[iy+q 2 + i/_(fc - g) 2 ] 
v-q 2 [{v + q 2 +V-(k- q) 2 } 2 + B 2 (k - q) 2 

D-(k-q) 



V-Bl(k - qf 



finite, 



(43) 



since in the long wavelength limit B 2 (k — q) 2 dominates over {v+q 2 + v~ (k — q) 2 } 2 - We have considered only the real 
part of S_(fc) (only this part will renormalise v~). However this does not produce any infra-red singular correction 
to i*2 and hence we ignore it. Note that it is the presence of the kinematic waves which makes this diagrammatic 
contribution infrared finite and thus ignorable. The second diagram in Fig. (10) also has a similar finite form and is 
ignored again. 




FIG. 11. One-loop contribution to C2{k,ui). A line refers to a response function and a line with a small filled circle refers to 
a correlation function . 



In Fig, 11 the diagram has the form 



C-(k, U )~- 



dq q 2 (k~q) 2 



fc 4 J q 2 (k - q) 2 [v-q 2 + v+{k - q) 2 - iB Q (k - q)} ' 



(44) 



It is easy to see that this is not infrared divergent. So we ignore all corrections to G_(fc, to). We immediately obtain 
z_ = 2, X- — 1/2- Note that the presence of the wave term Bgk in the inverse propagator was again crucial to this 
analysis. We see from the foregoing analysis that 



GZ 1 {k,w) ~iu + k 2 , C-(k,u = 0) 



/fc 4 ' 



(45) 



We now calculate the exponents of h + . As stated earlier we ignore corrections to A2 and A3 as they are all finite. 
There are however, diverging one-loop corrections to both C+ and E + which are shown in Fig.OI and Fig.|l3|. 





(a) 0» 
FIG. 12. One-loop contribution to C+(k,u)). A line refers to a response function and a line with a small filled circle refers to 

a correlation function. 



1G 



(a) (b) 
FIG. 13. One-loop contribution to E+(fe, w). A line refers to a response function and a line with a small filled circle refers to 

a correlation function. 

In both Figs. |l2| and [l^ diagrams (a) comes from the KPZ nonlinearity. If these diagrams were the only diagrams, 
one would have obtained KPZ exponents:z + = 3/2, \+ = 1/2- We however notice that the diagrams (b) in Figs.[l2 



and 13 are as strongly infrared (IR) divergent as diagrams (a). Let us examine diagram (b) in Fig.h3 in detail. The 



integral is given by (using the fact that \- = 1/2, z_ = 2) 



q 2 (k-q) f (k-q)[v-q 2 + v_(k-qf 



dQdq (n 2 + l£g*)[-*n + iB a k + V-(k- q) 2 ] ~ k J dq ^[B 2 k 2 + {v-q 2 + i/_(Jfc - q) 2 } 2 } 

*-WXmi ~ - / % + °( fc2 /, 4 ) + ~ ^ (46) 
v- \B z k l + Ais±q A \ j/_ qi 

Note that the presence of the kinematic wave was crucial again in this evaluation: The integral in (^) diverges as 
J dq/q 2 if the external wavenumber k is set to zero. At nonzero k, the integral is controlled by the presence of the 
wave term B 2 k 2 in the denominator, and scales as fc -1 / 2 , not as k~ x as might have naively been expected. Thus 
z + = 3/2 is unaltered by the second diagram. Similarly we consider diagram (b) of Fig.[l3] (using z + = 3/2): 

1 f 1 f 1 1 

1 ~ P J dndq (W + i/2 g4)[(n _ Bo k) 2 + v 2 _q^\ ~ iu_k?B k J dm ^B k-iv_q 2 ~ B k-iv_q^ ~ k ' (4?) 

which is as divergent as the diagram (a) in Fig.|l3|. Thus \+ = 1/2 also remains unaltered. It is easy to see that with 
these self-consistent response and correlation functions, one-loop corrections to A2, A3 do not diverge. Hence ignoring 
A2, A3 from our selfconsistent calculation is justified. 

So far we have worked in the comoving frame of i.e, with Eqs. (j33"[). Let us now go to the comoving frame of 
h + : We work with Eqs. (|3^). The relevant diagrams are same as given in Figs, [l^ and [l^. It is easy to carry out 
a self-consistent analysis again on Eqs. ({j2|). The only difference is that now C_|_(fc,aj) is peaked at oj = whereas 
C_(fc, oS) is peaked at lu = —B a k. Here again there is no singular correction to the reponse and the correlation 
functions of the field Thus z_ = 2 and x- = 1/2- There are however diverging corrections to v\ and D± which 
are identical to those obtained in our calculations in the comoving frame of Reassuringly, we find again same 
exponents, i.e. z + = 3/2, x+ — 1/2, as we should. Thus the effective equations in the left and right going frames are 



dh+{k) 
dt 



iB kh + = v + k 3/2 h + (k, t) + .... 
dh^{k) 



Of 



v_k 2 h-(k,t) + .... (48) 



and 



dh^{k) 
dt 



iB kh_ = v^k 2 h-(k,t) + .... (49) 



The renormalized equations ( pq ) and ( |4S| ) arc thus connected by a Galilean transformation like the bare equations 
(p3|), and (p2), and the exponents are frame- independent. 



C. Renormalization group analysis 



In Section VB we have shown how a self-consistent mode coupling treatment of our model for certain parameter 
values leads to weak dynamical scaling, i.e., distinct dynamical exponents for the two fields. In this Section we recast 
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these results within a perturbative dynamical renormalization-group (DRG) framework, and show how the difficulties 
posed for the DRG by weak dynamical scaling can be overcome. 



The diagrams are same as shown in Section VB. We work with Eqs. (|32j). As usual, the bare diagrams for 
z/_, £)_, A2 and A3 do not diverge. The diagram (a) in Figjl^ diverges as ~ 2±^k 2 J % ~ -±^fc 2 x i, whereas 
the diagram (b) in Fig.O diverges as \2\3D-k 2 J .„ , ^| ? ~ A 2 A 3 Z)_fc 2 x -4=. Thus diagram (b) can be ignored 

^ o "T" I/ — *? V fe 

in comparison with (a) in the small ft limit. Similarly, diagram (b) in Fig.O is less divergent than the diagram (a). 
Thus all relevant diagrams generating dominant singular corrections to the bare response and correlation functions of 
the field h+ are of KPZ type (with identical symmetry factors). 

Starting with a cutoff wavenumber A, eliminating modes with wavenumbers between Ae~ si and A, rescaling so that 
the cutoff wavenumber is once more A, and passing to the limit 51 — > 0, we obtain the differential recursion relations 

dv + r o , #1 

*B±=D + [z-\-2 X+ + ^ (50) 

and 

^ = D4z-l-2 X -}, (51) 

where v+, D + and £>_ are now functions of i, and g = D + \\/v^_ is the dimensionless coupling constant. By using 
g = 2 at the stable RG fixed point |l0| , we obtain 

^jt = !/+[*- 3/2], (52a) 

= „_[*- 2], (52b) 

= £>+[«- 1/2 -2 X +], (52c) 

^=U_[z- l-2x-]. (52d) 

Now given the exponents z, x+,X— > the equations of motion ( p2[ ) are supposed to be invariant under the scale 
transformations x — > 6x, t — > 6 z i, /i+ — * b x +h+, h- — ► b x -fi-. From Eqs. (|50| ) we get 2 = 3/2, x+ = 1/2 whereas 
Eqs. ( |5l|) give z = 2,%_ = 1/2. Which value of z should we choose? The choice of z = 3/2 keeps i/i fixed under 
rescaling, but v^{l) ~ e~'/ 2 as the RG fixed point is approached (remember e' is like a length scale). On the other 
hand if we choose z = 2 then ^ 2 is fixed under rescaling but v\ ~ e'/ 2 as the fixed point is approached. Moreover, 
z = 2 implies X- = 1/2 giving C_(a;,0) ~ x in agreement with the previously obtained results, and \+ = 1/4 
suggesting C + (x,0) ~ a; 1 / 2 , which, of course, is wrong. Similarly, the choice z = 3/2 gives correct spatial dependence 
for C + (x,0) ~ a;, but gives an incorrect result for C_(x, 0). 

To resolve this difficulty, let us first study the seemingly trivial case of two totally decoupled fields: 

d<t>i d 2 (j>! 

-dT = fll ^ + m ' (53a) 

Q ^ 2 92 ^ . ftOU 

with (??i(0, 0)rjj(x, t)) = 2^/iSijS(x)d(0) for i = 1,2 and no sum on repeated indices. These are linear equations and 
hence exponents can be found exactly. In particular we know 

(<MM)^i(0,0))~z/iOr7i), (54a) 
(0 2 (M)(MO,O)) ~x 3 / 2 (z 4 A)> (54b) 

i.e., we have Xi = 1/2, X2 = 3/2, z\ = 2, z 2 = 4. There are obviously no diagramatic corrections to any of the 
parameters. Here the flow equations are 
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^ - -2]: < r >- r >"'> 

'^f - I - 2\,]: 



= // 1 \z — Z 


• 

1 ! 


= 7! [z - 1 


-2xi 


= M2[z - 4 


!; 


= 72 - 1 


-2X2 



(•>•><■) 

^=72^-1-2x2]. (55d) 

A choice of z = 2 and xi = 1/2 keeps /^i and 71 fixed. This immediately tells us that equal-time correlation function of 
the field cf>\ scales as Ax 2xi . The coefficient A, a function of 71 and /ii is scale independent. However a naive use of this 
value of z (i.e., z = 2) leads to wrong conclusion about the spatial dependence of the equal time correlation function 
of the field of <\>2- z — 2 gives xi = 1/2 suggesting (02(0, 0)02 (ir, 0)) ^ x which is wrong (the correct dependence 
is a; 3 which is known exactly) . The reason is that the choice z = 2 makes the solutions of the flow equations for \i2 
scale-dependent: fi2 ~ e~ 21 ~ k 2 . However, one can still extract the correct behaviour of the correlation function of 
4>2 in the above example if one accounts for the scale dependence /12 explicitly: An explicit construction gives 

(Mk,")<h(-k,u))= * u + . ] ,J , (56) 

H2K X IUJ + ll2K z —IUJ + H2K 

which gives (by using scale dependent ^2) (<^i(0, 0)<j>2(x, 0)) ~ x 3 which is the correct answer! This suggests that to 
make sense of out the RG flow equations in presence of different dynamic scalings in a coupled system, one has to take 
care of the scale dependent diffusion coefficients while calculating the spatio-temporal behaviour of the correlation 
functions correctly. Let us review our RG results obtained from the flow equations ( |50|) and (|5l| ) in view of our 
previous analysis: We choose z = 2, X- = 1/2 which gives x+ = 1/4 and makes v\ ~ fc -1 / 2 . Using this scale 
dependent coefficient we correctly obtain 

C+(x,t) =xf-[x 3 / 2 /t}. (57) 

This suggests that an RG treatment, suitably modified, can be applied successfully to a problem of weak dynamic 
scaling provided scale dependent dissipation coefficients are taken into account while constructing the correlation 
functions. 



D. Strong dynamical scaling 

In Section H we have seen that for model equations |8|both the dynamic exponents turn out to be 3/2 (i.e., KPZ- 
like). This is easy to understand analytically: In this case both the equations have KPZ-like nonlinearities in addition 
to the non-KPZ ones. However, it is easy to see that for both the fields, due to the presence of the waves, in any 
frame, non-KPZ diagrams are less singular than the corresponsing KPZ-like diagrams. Thus the dominant singular 
corrections to the response and correlation functions for both the fields are KPZ-like, making x+ = X- = 1/2 and 
z + = z_ = 3/2. 



VI. SUMMARY 



This paper is a study of the nature of spatiotemporal correlations in a coupled-field driven diffusive model (the 
LR model [|]]^]), in the phase in which it displays kinematic waves. Our results include a demonstration of pairwise 
balance for certain parameter values and, hence, a proof that the steady state has a product measure in that range 
of parameters. Most important, we have been able to show that the fields corresponding to the two eigenmodes 
of the linearized version of the model were characterized by two different dynamic exponents, although the fields 
themselves are (nonlinearly) coupled. This is the first demonstration of such weak dynamical scaling in a model with 
two fully coupled fields (as distinct from the model of |3|,|l8],^0| , where such an effective decoupling was found only in 
semi- autonomously coupled systems in which one of the fields evolve independently). We have been able to show this 
numerically, through Monte Carlo simulations on a lattice model, and analytically, using self-consistent perturbative 
calculations as well as symmetry arguments, in the corresponding continuum stochastic PDEs. Outisde this special 
subspace of parameter space, the model exhibits normal, strong dynamic scaling. We also discuss and largely resolve 
the technical difficulties in applying the dynamical renormalization group when weak dynamic scaling prevails. 
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